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ABSTRACT 

Out of the known transiting extrasolar planets, the majority are gas giants orbiting 
their host star at close proximity. Both theoretical and observational studies sup- 
port the hypothesis that such bodies emit significant amounts of flux relative to the 
host star, increasing towards infrared wavelengths. For the dayside of the exoplanet, 
this phenomenon typically permits detectable secondary eclipses at such wavelengths, 
which may be used to infer atmospheric composition. In this paper, we explore the ef- 
fects of emission from the nightside of the exoplanet on the primary transit lightcurve, 
which is essentially a self-blend. Allowing for nightside emission, an exoplanet's tran- 
sit depth is no longer exclusively a function of the ratio-of-radii. The nightside of an 
exoplanet is emitting flux and the contrast to the star's emission is of the order of 
~ 1CP 3 for hot-Jupiters. Consequently, we show that the transit depth in the mid- 
infrared will be attenuated due to flux contribution from the nightside emission by 
~ 10~ 4 . We show how this effect can be compensated for in the case where exoplanet 
phase curves have been measured, in particular for HD 189733b. For other systems, 
it may be possible to make a first-order correction by using temperature estimates of 
the planet. Unless the effect is accounted for, transmission spectra will also be pol- 
luted by nightside emission and we estimate that a Spitzer broadband spectrum on 
a bright target is altered at the 1-a level. Using archived Spitzer measurements, we 
show that the effect respectively increases the 8.0/im and 24.0/im transit depths by 
1-a and 0.5-er per transit for HD 189733b. Consequently, we estimate that this would 
be ~5-10cr effect for near-future JWST observations. 
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1 INTRODUCTION 

The field of exoplanet ary science was g alvanized by the dis- 
covery of t he first trans it ing p lanet by ICharbonneau et al.l 
||2000T) and iHenrv et all (|2000h . The transit discovery was 
made in the visible wavelength range and for the subse- 
quent few years this was established as the normal pr ac- 



tice i n later observ a tions a nd surveys (e .g. Brown et al.l 



i|200ll h iBakos et all (|2004l h iPollacco et all (|2006l )h In re- 
sponse to the growing need for accurate parametriza- 
tion of lightcurves, several authors produ ced equations 
model li ng the lightcurve behavio r inclu di ng ISackett et al.l 
Jl999h. ISeager fc Mallen-Ornelas! (|2003h . iGimenezI (|2007f ) 
and lKippind (|2008h . 

In the last two years, the value of infrared measure- 
ments of transiting systems has become apparent with nu- 
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merous pioneering detectio ns; emission from a transiting 
planet (iDeming et al . 2005), emission from a non-transiting 
planet (Harrington et al.l 120061 ). an exoplanetary spectrum 
IjGrillmair et al. 120071). detection of w ater in the atmosphere 



of HP 189733b l|Tinetti et al.ll2007bl ). me thane ilSwain et all 
120081 ) and more recently carbon dioxide (I Swain et alj 20091 ) . 
More details on the use of transmission spectroscopy as 
a tool for detect i ng mo lecula r species can be f o und i n 
ISeager fc Sasselovl |2000l ) and iTinetti fc Beaulieul (|2009h . 
With JWST set to replace HST in the next decade, we 
can expect an abundance of high- precision infrared transits 
to be observed in order to dete ct more molecular species, 
perhaps including biosignatures l|Seager et al.ll2005l ). In this 
paper, we discuss the consequences of significant nightside 
planetary emission on precise infrared transit lightcurves. 
Nightside emission renders one of the original assumptions 
of transit theory invalid: 'The planet may be treated as a 
black disc occulting the flux-emitting star'. In the case of 
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hot- Jupiter systems, the nightside of the planet is now also 
flux-emitting. This additional flux can be considered as a 
blend, but from the planet itself, i.e. a self-blend. 

Conceptually, it is very easy to see that this will cause 
mid-infrared transit depth measurements (and to a lesser 
degree in the visible and near-infrared range) to become un- 
derestimates of the true depth. The reason is that there are 
two sources of flux, the star and the planet, and only one of 
these is being occulted, whilst the other is the blend source. 
Note that we define the true transit depth as a purely ge- 
ometric effect determined by the ratio-of-radii squared i.e. 
dgco = (Rp/R*) 2 = p 2 . In this work, we will derive expres- 
sions estimating the amplitude of the effect, propose meth- 
ods for correcting the lightcurves and apply them to two 
cases where the nightside emission of an exoplanet has been 
determined in the mid-infrared. 



2 DERIVATION 
2.1 Depth dilution 

Let us define the total out-of-transit flux surrounding a tran- 
sit event, as shown in figure 2a, to be given by: 



Fout 



F t + Fp, nig ht 



(1) 



where F* is the total flux received from the star and 
Jp,ni g ht is the total flux received from the nightside of the 
planet, over a time interval of St. Let us assume that the star 
is a uniform emitter and that both the stellar and planetary 
total flux are invariable over the timescale of the transit 
event. We may then write down the flux during a transit as 
a function of the ratio of the radii, p: 
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F* + F P . 



night 



(2) 



Note how the flux of the star has been attenuated as a 
result of the eclipse but the planetary flux is still present. 
The observed transit depth in the flux domain, d b s , is de- 
fined by: 

T F ut ,tra Firi ,tra 



Fout,tTPL 

F, 



(3) 

dobs = ( " * ) ■ rfgco (4) 

v r* + -TP, night / 

Where d gco is the geometric transit depth. In the case 
of Jp jn i g ht — > 0, we recover the standard equation for the 



depth being d a 



p . For cases where the nightside 



flux of the planet is non-negligible, the transit depth will 
therefore be affected. We may re- write equation (4) as d geo = 
Snight x dobs where we define: 



F,+F, 



night — 



P.night 



(5) 



The blend source is at a much cooler temperature than 
the host star and thus the contrast between the two bodies 
is greater at infrared wavelengths. Consequently, the self- 
blending has a significant effect on infrared measurements 
(e.g. see §5.1) but a much lower impact on visible wave- 
lengths (e.g. see §4.4). This makes the inclusion of such an 



effect paramount since visible and infrared measurements 
must be considered incommensurable unless this systematic 
is corrected for. 



2.2 The consequences for other parameters 

The largest effect of a blend is to underestimate the transit 
depth. However, we evaluate here the effect of the night- 
side blend, or indeed any kind of blend, on the other 
lightcurve parameters. For simplicity, we consider here 
a circular orbit and so we m ay use the expressions of 
ISeager fc Mallen-Ornelasl (|2003h . Assuming a blend factor 
given by B we have d &s = P 2 /B. Sea ger fc Mallen-Ornelasl 
(|2003T) derived expressions for retrieving the impact param- 
eter squared, b 2 , and the semi-major axis (in stellar radii) 
squared, (a/R,) 2 , as a function of tr, tF and d, where tr 
and t F are the l st -to-4 th and 2 nd -to-3 rd contact transit du- 
rations respectively. We may calculate the retrieved value of 
b 2 and (a/R*) 2 in the case where a blend exists but we have 
negated it in our analysis . Using equations (7) and (8) from 
ISeager fc Mallen-Ornelasl (|2003l ). we find: 



-"derived 
1.2 



p 2 - (1 + p 2 - b 2 )VB + B 
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bje^ed = b 2 + 1(1 - b 2 - p z )(B - 1) + 0[(B - I) 2 ] (6) 



(B-l) 



( a /R*)derived — {a/R*) 2 

(l-b 2 -p 2 )((a/R,) 2 -(l+p) 2 ) 
2((l+p) 2 -& 2 ) 
+ 0[(B-1) 2 ] (7) 

Equations (6) and (7) imply that negating the blend- 
ing factor causes us to overestimate the impact parameter 
and underestimate a/R,. It is important to recall that the 
lightcurve derived stellar density is found by taking the cube 
of a/R* and therefore will exacerbate any errors at this 
stage. We note that both equations give the expected re- 
sults for B = 1, i.e. no blend source present. 

These expressions have been derived assuming no limb 
darkening is present, which is typically a very good approx- 
imation for the wavelength range we are interested in. How- 
ever, in reality the incorporation of limb darkening is easily 
implemented and demonstrated later in §5. Therefore, equa- 
tions (6) and (7) should not be used to attempt to correct 
parameters derived from fits not accounting for nightside 
pollution, rather they offer an approximate quantification 
of the direction and magnitude of the expected errors. 



3 COMPENSATING FOR THE EFFECT 
3.1 Empirical method 

In the previous section we saw how the geometric transit 
depth, dgoo, and the observed transit depth, d b s , are re- 
lated by the factor B n ight- Fortunately, B n ight is an observ- 
able and can be obtained through measuring the phase curve 
of an extrasola r plan et (for the first measured example, see 
iKnutson et al.l (|2007l )). With such a measurement, the dif- 
ference between the day and nightside fluxes may be deter- 
mined and thus Fp in ight/F* can be calculated. 
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Another possible method is to measure the secondary 
and primary transits without the intermediate phase curve 
information, which would require instruments with ex- 
tremely stable calibration. For a very inactive star, a highly 
calibrated instrument could, in principle, measure the night- 
side flux by measuring just the primary and secondary 
eclipses. This would equate to an absolute calibration ac- 
curate to a fraction of the difference between the day and 
nightside fluxes, estimated to be ~ 10 -3 for HD 189733b, 
over one half of the orbital period. We therefore estimate 
calibration requirements to be at least ~ 10~ 4 between a 
~ 30 hour period. Whilst this would be a moot point for 
staring telescopes like Kepler and CoRoT, it is the infrared 
telescopes of Spitzer and JWST that are most heavily af- 
fected by nightside pollution and these telescopes frequently 
slew around looking at different patches of the sky. After the 
slewing we require the target to be at the same centroid posi- 
tion to within a fraction of a pixel. Whether the whole phase 
curve or simply the eclipse-only observations are made, the 
same method may be used to correct both for the effects of 
nightside pollution. 

A phase curve time series is typically normalized to F* , 
as shown in figure 1, which is in contrast to a normal tran- 
sit measurement, which is normalized to F* + Fp jn i g ht i.e. 
the local out-of-transit baseline (see figure 2a for illustra- 
tion). Therefore, for a phase curve, the stellar normalized 
flux immediately before and after the transit event is equal 
to Bnight, as shown in figure 2b. ft is possible that shifted 
hot spots on the planetary surface could cause an inequality 
between the pre and post transit baselines, but in practice 
the net effect of nightside pollution is very well accounted 
for by averaging over this time range. 

In order to correct a primary transit lightcurve, we need 
to modify the normalization. In figure 2, we show the two- 
step transformation which can achieve this. We consider ini- 
tially normalizing a lightcurve using the local baseline as 
usual for such measurements, as shown in figure 2a. After 
this, the two-step correction may be performed, provided the 
observer has knowledge of B n i g ht- The whole process may be 
summarized by the following (also summarized in figure 2): 

1. Normalize fluxes to local out-of-transit baseline, as 
usual. 

2. Multiply all flux values by Bnight- 

3. Subtract (B n ight — 1) from all data points. 

In practice, these steps are incorporated into the 
lightcurve fitting algorithm directly. In the case of using 
Monte Carlo based routines for error estimation, £> n i g ht may 
be floated around its best-fit value and corresponding uncer- 
tainty. An example of this method is shown in §5.1 and §5.2 
for the planet HD 189733b. Defining ij, U ncorr as the locally 
normalized flux measurement of the j th data point, we may 
explicitly write down the corrected data point as: 



presence of ellipsoidal variations of the star, which would 
mix t he phase curve signature. For example, Wel sh et al.l 
(2010) detected ellipsoidal variations in HAT-P-7. Such sig- 
nals peak at orbital phases of 0.25 and 0.75, whereas a phase 
curve should peak close to orbital phase 0.5, but can be off - 
set by a small factor due to hot spots. IWelsh et al.l (|2010h 
provide a detailed discussion of modelling both signals and 
although ellipsoidal variations complicate the analysis, they 
certainly do not undermine it. 

3.2 Semi-empirical method 

To correct for the effect of nightside pollution, in an accu- 
rate way, we have proposed using phase curve information to 
obtain Bnight , which requires many hours of telescope time. 
Further, the phase curve should be obtained at every wave- 
length simultaneously and for every epoclQ one wishes to 
measure the transit event at, in order to be sure of a com- 
pletely reliable correction. We label this resource-intensive 
method of correcting for the effect as the 'empirical method'. 
However, we appreciate that obtaining phase curves at every 
wavelength for each epoch is somewhat unrealistic and pro- 
pose a 'semi-empirical method' of achieving the same goal 
with far fewer resources. 

We propose that observations of the phase curve are 
made at several wavelengths in the infrared; for example 
the IRAC and MIPS wavelengths of Spitzer are very suit- 
able but most of these channels are unfortunately no longer 
available. With only a few measurements of the nightside 
flux, a prudent approach would be to assume the star and 
planet behave as blackbody emitters as a first approximation 
and extrapolate the emission to other wavelengths which 
are missing phase curve observations (an example of this 
will be provided in §7). Should higher resolution spectra of 
the nightside emission become available in the future, we 
may construct of a more sophisticated model as appropri- 
ate. Transit observations at different wavelengths may then 
interpolate/extrapolate the model template to estimate the 
magnitude of the effect and apply the required correction 
at any wavelength. This allows us to estimate the blending 
factor at all wavelengths and times. 

Long-term monitoring of the planet may also be neces- 
sary in order to ascertain the presence or absence of tempo- 
ral variability in the system. Ideally, this may be achieved 
by obtaining phase curves of the exoplanet to measure the 
nightside flux at regular times. More practically, it could 
be done by measuring only the secondary eclipse, which is 
the dayside of th e exoplanet, at regular times (for example 
lAgol et"al] (|2009T l). Any large changes in the nightside flux 
are likely to be correlated to large changes in the dayside 
flux too, assuming a constant energy budget for the planet. 
This second approach would reduce the demands on tele- 
scope time by an order of magnitude or more. 



1 J .corr 

= B night J- j .uncorr - (Bnight " 1) (8) 

We also briefly mention here that the transformation 
operations on the lightcurve time series will not only change 
the transit depth but also provides a more physical transit 
signal and thus we expect a slightly lower \ 2 m the final 
fitting, as indeed is seen later in tables 1 and 2. 

One caveat with the described method is the possible 



3.3 Non-empirical method 

The final method we propose here is the least accurate but 
requires the fewest resources to implement. If a planet has 

1 Each epoch should be done in case temporal variability exists 
in the system 
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Figure 1. Predicted phase curve of a hot-Jupiter with a hot dayside and cooler nightside. Note how all fluxes are normalized to the 
star-only flux. 



recently been discovered, no phase curves or even secondary 
eclipses may have been obtained yet. Concordantly, the only 
avenue available is to estimate the temperature of the night- 
side through either a simple analytic estimation or a dynam- 
ical model of the atmosphere, although the latter may be ex- 
cessive given the absence of any observational constraints. 
We illustrate here how the temperature may be quickly es- 
timated in such a case. 

The nightside temperature may be estimated by assum- 
ing the dominant source of heating is from the incident stel- 
lar flux. In this case, the only unknown factors affecting the 
nightside temperature are the re-distribution of heat factor, 
/ and the Bond albedo of the planet, Ab- First-order esti- 
mations of these values can be made based upon empirical 
upper limits and measurements of other planets and atmo- 
spheric models. The following expression may be used as a 
first-order estimate of the planet's brightness temperature: 



P,hcmisphere 



(A) 



+ 



L ln 



AttR p ob£ 



1/4 



(9) 

Where T* is the effective temperature of the host star, 
7?* is the stellar radius, ap is the orbital semi- major axis of 
the planet, / is the distribution of energy to the hemisphere 
in question, L- m % is the luminosity of the planet from inter- 
nal heat generation (e.g. tidal heating, radioactivity), <tb is 



the Stefan-Boltzmann constant and Rp is the radius of the 
planet. 

For this calculation, Lint is generally assumed to be zero 
unless large tidal forces are expected as a result of a highly 
eccentric orbit, for example. All of the other quantities are 
typically measured except for / and Ab- Choices for these 
values may come from atmospheric models or experience 
with other exoplanets. 



4 COMPARISON TO OTHER EFFECTS 
4.1 Starspots 

Starspots have b een observed with i n the transit events i n 
several cases; e.g. iPont et all i|2007l ). lDittman et al.l l|2009l) . 
They typically have been observed to have a radius of a few 
Earth radii and are estimated to have temperatures from 
100K to 1000K cooler than the rest of the stellar surface. 
When a planet passes over a starspot, it results in an in- 
crease in relative flux within the transit signal which is eas- 
ily identified. If one assumes that the only starspot is the 
starspot which has been crossed, then accounting for the ef- 
fect is quite trivial and may be incorporated in the lightcurve 
modelling. 

What is much more troublesome are out-of-transit 
starspots for which we have no direct evidence. The pres- 
ence of out-of-transit starspots will cause the transi t depth 
to appear larger, in general. As a typical example, ICzeslal 
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Figure 2. Illustration of the three-stages involved in the corrective procedure, to compensate for the effects of nightside pollution. 
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(2009) estimate that the effect can cause underestimations 
of the planetary radius by a fraction of ~ 3% for CoRoT-2b, 
which is a 1.6% change in the transit depth. This effect is 
larger than the nightside pollution effect by a factor of 5-10. 

However, these effects can only be present for spotty 
stars whereas nightside pollution simply requires a hot 
planet. Also, the spot coverage of a stellar surface varies pe- 
riodically giving rise to an ultimately regular pattern which 
may therefore be corrected for. In contrast, the nightside 
pollution effect is not periodic, it is a constant offset in a 
single direction. Further, prior information such as a phase 
curve or a dayside eclipse places strong constraints on £? n ight 
and therefore the estimation of this parameter is not an is- 
sue. Therefore, even for spotty stars, no-one would propose 
ignoring the effects of blending induced by a nearby compan- 
ion star and so it can be seen that negating the self-blend 
of the planet would also be folly. 

4.2 Temporal variability of Fp and F* 

If the stellar flux or the nightside emission of the planet ex- 
periences temporal variations, then we would expect d b s to 
also change over time. We will here estimate the magnitude 
of this effect. The expected changes in nightside emission 
has not been studied in as much detail as that for the day- 
side, b ut we expect the magn itude of variations to be very 
similar. [Rauscher et al.l(|2007T ) have used shallow-layer circu- 
lation models to estimate variations in the dayside emission 
at the 1%-10% level. This is con sistent with the ob servations 
of HD 189733b's dayside flux bv lAgol et ail (|2009h who mea- 
sure variations in the dayside below 10%. Typical stellar flux 
changes are at the 1% level and so the ratio Fp/F* is more 
likely to vary due to the planet than the star. 

For the case of HD 189733b, we later show that Bnight = 
1.002571 at 8.0/xm. If F night increased by ±10%, this would 
correspond to 5 n i g ht = 1.002828 causing the transit depth 
to vary by 0.0006%. This would be around an order of mag- 
nitude below Spitzer's sensitivity but could be potentially 
close to a 1-a effect for JWST. Nevertheless, the effect is 
sufficiently small that it is unlikely to be significant in most 
cases. 



4.3 Limb darkening 

The nightside pollution effect is generally only relevant for 
hot-Jupiters at infrared wavelengths. As we move towards 
the 10-30 (jm wavelength range, the effects of limb darken- 
ing become negligible. The curvature of the transit trough 
is essentially flat. However, th e limb of the star w i ll pos - 
sess a more compli cated profile (|Orosz k, Hauschildtl |2000); 
Ijeffers et alj |2006)) and this could potentially introduce er- 
rors into the fitting procedure. It is generally prudent to in- 
clude even the very weak limb darkening effects when mod- 
elling such transits. 

As a result, equations (6) and (7) should not be used 
to attempt to correct lightcurves which were fitted without 
nightside pollution. They do, however, offer a useful approxi- 
mate quantification as to the direction and magnitude of any 
errors. A comparison between the predictions of (6) and (7) 
and the exact limb-darkened nightside-polluted lightcurve 
fits is given later in §5.1 for the example of HD 189733. 



A discussion of this method required to produce this exact 
modelling is given in §3.1. 

4.4 Significance at visible wavelengths 

We briefly consider the value of including the nightside 
pollution effect at visible wavelengths, in particular for 
the Kepler Mission. iBorucki et al.l l|2009h recently reported 
visible- wavelength photometry for HAT-P-7b which exhibits 
a combination of ellipsoidal variations and a phase curve 
l| Welsh et alJl201Ch . as well as a secondary eclipse of depth 
(130 ± ll)ppm. HAT-P-7b is one of the very hottest tran- 
siting cxoplanets discovered, so it offers a useful upper-limit 
example. For the purposes of nightside pollution, the maxi- 
mum possible effect would occur if eclipse was both due to 
thermal emission alone and efficient day-night circulation. 
In this hypothetical scenario which maximizes the night- 
side pollution, we would have B n i g ht = 1.00013. In the 
case of HAT-P-7b, the t ransit depth was reported to be 
dobs = (6056±47)ppm bv lWelsh etafl l|2010l ), implying that 
the geometric transit depth is larger by 0.79ppm, or 0.017-ct. 
Therefore, as expected, visible wavelength transits will not, 
in general, be significantly affected by nightside pollution for 
even Kepler photometry. 



5 APPLIED EXAMPLE - HD 189733B 

5.1 Spitzer IRAC 8.0/im measurement 

We will here provide an example of the empirical method 
of compensating for nightside pollution. We used the cor- 
rec ted data of HD 18973 3's phase curve at 8.0/im, as taken 
by iKnutson et al.l (|2007l ) (obtained through personal com- 
munication). We applied a median-stack smoothing function 
to the lightcurve with a one-minute window in order to iden- 
tify the eclipse contact points. We find the flux of the star by 
taking the mean of fluxes between the 2 nd and 3 rd contact 
points during the secondary eclipse, weighting each point 
by the reported error. The standard deviation within this 
region is divided by the square root of the number of data 
points to give us the error on the mean. All fluxes are then 
divided by the derived stellar flux and the error on each flux 
stamp is propagated through, incorporating the error on the 
stellar flux estimate. 

In order to determine the nightside flux, which is not 
the same as the minimum flux, we adopt a baseline defined 
as 30 minutes before 1 st contact and 30 minutes after 4 th 
contact and find a mean of B night = 1.002571 ± 0.000048. 
The average rms in this baseline is 0.65 mmag/minute. If 
it were possible for the nightside of the planet to induce a 
secondary eclipse, as the dayside does, we would measure 
a secon dary eclipse nightsid e dept h of (0.256 ± 0.023)%, 
whereas ICharbonneau et al.l (|2008l ) report a dayside sec- 
ondary eclipse depth of (0.391 ± 0.022)%. 

As discussed in §3.1, ellipsoidal variations can also be 
responsible for out-of-transit flux variations and can be po- 
tentially confused wit h the phase curve . For HD 189733, 
we use equation (1) of IPfahl et al.1 (|2008l ) to estimate an el- 
lisoidal variation amplitude of 2.2ppm. Given that the phase 
curve exhibits a variation 1350ppm amplitude, ellipsoidal 
variations can be neglected for the rest of this analysis. 
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We now produce two fits of the lightcurve: 1) no blend- 
ing factor 2) blending factor £> n i g ht included. Each lightcurve 
is fitted independently assuming a fixed period of P = 
2.2f85733 days and zero orbital eccentricity. The results of 
the fits are displayed in table f. 

For the fitting, we use a Markov Chain Monte Carlo 
(MCMC ) algo rithm which employs the geo metric model of 
iKippind (|2008h . the limb darkening model of lMandel fc Agoll 
(2002) an d utilizes the li ghtcurve fitting parameter set de- 



scribed in lKippind (l20ld ): t c ,b 2 , T/R,, p 2 and OOT. We 
use f25,000 trials with the first 25 ,000 dis c arded for burn-in. 
Employing the code of I. Ribas, a lKuruczl (|2006l ) style atmo- 
sphere is used to interpolate t he four non-linear limb dark- 
ening co efficients JClard feoOO) , following the same method- 
ology of iBeaulieu et all l|20ldh , giving us d = 0.2790207, 
c 2 = -0.1506885, c 3 = 0.0779481 and c 4 = -0.0087653. We 
use the same local baseline as defined earlier, constituting 
22382 data points and assume a circular orbit. At this stage, 
no outliers have yet been rejected but we proceed to fit the 
unbinned lightcurve. We take the best-fit lightcurve and sub- 
tract it from the data to obtain the residues and then look 
for outli er points. W e use the median-absolute-deviation 
(MAD) jGaussI (|l816l )) to provide a robust estimate of the 
standard deviation of the data, as this parameter is highly 
resistant to outliers, and find MAD = 2.92062 x 10~ 3 . Since 
there are 22382 points, then the maximum expectant depar- 
ture from a normal distribution is 4.08 standard deviations. 
Any points above this level are rejected, where our definition 
of the standard deviation comes from the MAD value mul- 
tiplied by 1.4826, as appropriate for a normal distributior0. 
This procedure rejects any points with a residual deviation 
greater or equal to 0.0176749, corresponding to 10 pointfl 
The new lightcurve is then refitted in the normal way 
and we present the best-fit value in table 1. We find the un- 
corrected lightcurve has a dept h of d b s = 2, 3824± .0061%, 
For comparison, we note that iKnutson et "all (|2009h report 
the fitted S.Ofim depth to be 2.387±0.006%, which is consis- 
tent with our value. Applying the correction due to night- 
side pollution, we find the geometric transit depth to be 
dg CO = 2.3884 ± 0.0061% meaning that the depth has in- 
creased by 0.006% corresponding to 1-er. From this example, 
it is clear that negating an effect which systematically biases 
transit depth measurements by ~ 1-a would be imprudent. 

Using equation (7), to first-order in (B — 1), we esti- 
mate that the impact parameter should be overestimated 
by 0.00069. Our fits reveal a very similar figure of 0.00060. 
Similarly for a/R* we predict an underestimation of 0.0034 
whereas the lightcurve fit finds 0.0062. As expected, the ef- 
fect of a blend is less pronounced on the other parameters. 

Based on the difference in collecting area, it is expected 
JWST will achieve a precision ~ 6.6 times greater than 
Spitzer, suggesting this nightside pollution effect will become 
significant at the ~ 5-10-sigma level for future infrared tran- 
sit observations of hot-Jupiters. Additionally, the binning of 
multiple Spitzer transits would raise the significance of the 



2 Although strictly a Poisson distribution, for 22382 data points, 
the distribution is very well approximated by a Gaussian 

3 We note that the data has already been cleaned of outlier mea- 
surements, which is why the MAD rejection criteria only identifies 
10 outliers from 22382 points 



effect. For example, lAgol et all (120091 ) reported seven 8.0/im 
transits of HD 189733b which, if globally fitted would inrease 
the significance of the nightside pollution effect to 2.6-cr. 
Such a large effect cannot be justifiably negated. 

5.2 Spitzer MIPS 24^m 

IKnutson et alj (|2009l ) measured the phase curve of HD 
189733b with the MIPS instrument onboard Spitzer about a 
year after the observations of the 8.0/im phase curve for the 
same system. Using the original normalized-to-stellar-flux 
unbinned data (personal correspondence with H. Knutson), 
we took the mean of data points ~ 1 hour either side of the 
transit event, which exhibit an rms of 2.1 mmag/minute. We 
combined to two baseline estimates to calculate a nightside 
rela tive flux of J3 ni ,h t = 1 .00438 ± 0.00025. 

IKnutson et al.1 (|2009h reported a 24/im transit depth of 
2.396 ± 0.027% and our own re-analysis of the data yields 
dobs = 2.398 ± 0.019%, where the fit has been performed us- 
ing the same methodology as for 8.0/xm, except we assume 
no limb darkening at 24/xm. As before, we apply the cor- 
rection due to nightside pollution and estimate a new 24/im 
transit depth of d gco = 2.409 ± 0.020%, which increases the 
depth by ~ 0.5-ct. Despite the absolute effect being larger 
than that at 8.0^im, the difference is fewer standard devi- 
ations away due to the much poorer signal-to-noise of the 
transit event itself. 



6 POLLUTION OF THE TRANSMISSION 
SPECTRUM 

In the standard theory of transmission spectroscopy, plane- 
tary nightside emission is assumed to be negligible and thus 
disreg arded ( e .g. as explicitly stated in the foundational the- 
ory of lBrownl (|200lh ^>. However, we have shown here that the 
effect noticeably changes the transit depth for high qual- 
ity photometry. Essentially, we posit that the 'traditional' 
transmission spectrum is in fact a combination of transmis- 
sion through the terminator and the self-blending caused by 
emission from the nightside. 

One subtle point is that the nightside pollution effect is 
not something which can be accounted for in the modelling 
of the transmission spectrum. It is generally useful to think 
of the nightside pollution effect as an astrophysical blend 
which happens to be related to the planetary properties. A 
transmission spectrum is typically found by fitting a transit 
lightcurve at multiple wavelengths and then fitting a spec- 
trum through the retrieved transit depths which models the 
planetary atmosphere. These routines usually make use of 
radiative transfer, chemical equilibrium, molecular line lists, 
etc to estimate the opacity of the atmosphere at each wave- 
length. However, attempting to increase the sophistication 
of these models would not accurately account for the self- 
blending scenario. Recall that each transit depth is obtained 
by fitting an eclipse model through the lightcurve time se- 
ries. Critically, it is at this stage where blending needs to be 
accounted for. The transit signal plus blend should be mod- 
elled as such from the outset due to subtle and quite intricate 
inter-dependencies between b, a/7?*, p 2 and the limb dark- 
ening. Thus it can be seen that attempting to incorporate 
the effects later on is far more challenging and completely 
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Table 1. Best-fit transit parameters for the HD 189733b 8.0fim primary transit lightcurve; data obtained bv \Knutson et a~l ] \200j) . Fits 
performed for the case of 1) typical normalization the local baseline 2) correction for the effects of nightside pollution. The number of 
data points is 22372. 



Method 


Depth, p 2 ,% 


T/s 


b 


a/R, 


i/° 


X 2 


(1) Local baseline 


2.3824 


5127.45 


0.66264 


8.9121 


85.7360 


22412.4691 


(2) Nightside correction 


2.3884 


5127.55 


0.66204 


8.9183 


85.7428 


22412.4605 


Uncertainity 


0.0061 


8.1 


0.0061 


0.050 


0.054 




(2) - (1) 


+0.0060 


+0.10 


-0.00060 


+0.0062 


+0.0068 


-0.0086 



Table 2. Best-fit transit parameters of the HD 189733b 24-Ofim primary transit lightcurve; data obtained bu \Knutson et a~l ] \2009\) . Fits 
performed for the case of 1) typical normalization the local baseline 2) correction for the effects of nightside pollution. Number of data 
points is 1198. 



Method 


Depth, p 2 ,% 


T/s 


b 


a/R, 


i/° 


x 2 


(1) Local baseline 


2.3980 


5072.527 


0.61425 


9.492 


86.290 


1260.9125 


(2) Nightside correction 


2.4085 


5072.505 


0.6131 


9.502 


86.300 


1260.9092 


Uncertainity 


0.019 


19.7 


0.003 


0.033 


0.022 




(2) - (1) 


+0.011 


-0.022 


-0.0012 


+0.010 


+0.010 


-0.0033 



unneccessary than simply fitting each transit lightcurve with 
a physically accurate model in the first place. 

In this section, we will estimate how different an exo- 
planet's transmission spectrum would appear with and with- 
out nightside pollution. In order to evaluate the magnitude 
of the effect, we will here consider a planet of similar type 
to HD 189733b. It is important to stress that the effects of 
nightside pollution will vary from case to case and the ex- 
ample we give here is indeed just for one example which is 
somewhat typical for an observed hot- Jupiter. Therefore the 
results here are only for a hypothetical, but typical, example. 

The real question we need to answer is how much does 
a planetary transmission spectra change due to nightside 
self-blending? We therefore need to generate two versions of 
the planetary transmission spectra, one including (figure 3a) 
and one excluding the effects of nightside pollution (figure 
3c), and then take the difference between the two (figure 
3d). 

For a description of the mo dels used to generat e the 
spectra, deta il s may be found in iTinetti et al.l (|2007aT ) and 
Tinetti et all (l2007bh for the transmission spectrum and 
Swain et air (2009) For the emission. Planet and star prop- 
erties are set to be that of the HD 189733 system. The 
model contains water, carbon dioxide and methane to pro- 
vide us with the effects of molecular species on nightside 
pollution. No carbon monoxide or hazes/clouds are included 
in our example. We note that the transmission and emis- 
sion models are good fits to the current available spec- 
troscopy /photomet ry data of HD 189733b in the NIR/MIR 
(|Swain et al.l 12009] ). The effects of w ater absorption ar e 
quantified with the BT2 water line list (Ba rber et al.| [2006). 
which characterises water absorption at the range of tem- 
peratures probed in HD 189733b. Methan e was simulated 
by u sing a combination of HITRAN 2008 (iRothman et al.l 
2005) and PNNL data-lists. Carbon dioxide absorption coef- 
ficients were estimated with a combination of HITEMP and 



CDSD-1000 (Carbon Dioxide Spectroscopic Databank ver- 
sion for high temperature applications; Tashkun & Perevalov 
(2008)). The c ontinuum was comp uted using H2 — H2 ab- 
sorption data l|Borvsow et alJfeoOll ). 

Generated spectra are always plotted in terms of the 
quantities determined with the lowest measurement uncer- 
tainty, namely (Rp/R*) 2 and (Fp/F t ), for the primary and 
secondary eclipses respectively. Transmission spectra which 
are plotted in units of Rp will cause the measurement uncer- 
tainties to be much larger since the error on R, must neces- 
sarily be propagated in such a recipe. In fact, the measure- 
ment uncertainties on a spectra plotted in units of Rp will be 
dominated by the error on R, since this property is typically 
determined to much lower precision. Consequently, statisti- 
cally significant molecular features would be overwhelmed 
by the artificially large error bars. 

Using the model described above, we first compute the 
transit depth from transmission absorption effects only (i.e. 
excluding nightside pollution) as visible in figure 3a. We then 
generate the dayside emission spectra for the same planet 
and make the assumption that the dayside and nightside 
emission spectra are identical (figure 3b). This assumption 
is unlikely to be true for the exact case of HD 189733b and 
really constitutes an upper limit, but we again stress that 
we are here only considering a planet similar to that of HD 
189733b and thus we are free to make this assumption for 
our hypothetical example. 

We combine the nightside emission and transmission 
spectra to produce a transmission spectra which includes 
the effects of the nightside, as seen in figure 3c. We then 
take the difference between the corrected spectrum and the 
one which excludes nightside emission. The resultant resid- 
ual spectrum is plotted in figure 3d. The residual spectrum 
reveals nightside emission affects the transmission spectra 
at the level of 6 x 10 -5 above 10 microns and very closely 
matches the behavior of the emission features, as expected. 
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Figure 3. Top left: A- Transmission spectrum of a hypothetical exoplanet similar to HD 189733b, generated considering the transmission 
through the terminator only. Top right: B- Emission spectrum from the nightside of the hypothetical planet. Bottom left: C- Transmission 
spectrum of the planet incorporating the pollution of the nightside emission. Bottom right: D- Residual between two transmission spectra. 
We conclude that not accounting for nightside emission would result in a 60-80ppm error in the transit depth. 



As we saw earlier, the magnitude of the effect is equal to 
the typical measurement uncertainty for a target like HD 
189733b with Spitzer. This supports our hypothesis that the 
nightside pollution effect is a ~1-<t effect for 8.0/^m Spitzer 
photometry. 



7 EXTRAPOLATING THE NIGHTSIDE 
CORRECTION FOR HD 189733B 

Only two measurements exist for the nightside flux for HD 
189733b (or indeed any other exoplanet) at 8.0/im and 
24^im, but several other primar y transit lightcurves exist 
in the mid-infrared wavelengths. iBeaulieu et al.l (|2008l) pre- 
sented 3.6/im and 5.8/xm measurements and Desert et al.l 
(2009) obtained photometry at 4.5/xm and 8.0/im. We will 
here estimate the nightside effect on the 3.6/im, 4.5/im and 
5.8/im channels. Currently, only two data points exist and 
this does not warrant us modelling the nightside emission 
spectrum in any more complexity than that of a blackbody, 
as a first order approximation. 

We first assum e that the star is b lackbody emitter with 
T» = 5040 ± 50K (|Torres et all 120091 ). At this point is ad- 
vantageous to consider only the emission per unit area from 
both the planet and the star, in order to avoid the effects of 



transmission through the planetary atmosphere. We there- 
fore define the flux per unit area of each object using: 



-7~P,night — -Fp,mght/(7T-Rp) 

T, = F*/(irRl) 

„ _ Fp, night _ Fp, night , 

'v-night — — r-, u>£ 



J 7 , 



(10) 

(11) 

(12) 



We then convert the B n i g ht measurements for S.Ofim 
and 24^im into 7^. n ight by using the two values of the cor- 
rected ratio-of-radii we showed in tables 1 and 2 and propa- 
gating the uncertainties. Using equations (6) & (7) and the 
relevant Spitzer bandpass response functions, we are able 
to estimate 7?. n ight for any given value of Tp !n i g ht. We allow 
this temperature to vary from 700K to 1700K in IK steps 
and numerically integrate the bandpasses to find x 2 a t each 
temperature, which we define as: 



+ 



3bs »F zaXc 

A(J- oba ) 



+ 



8.0Mm 
2 



A(j; bs ) 



(13) 
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Figure 4. Flux per unit area (i.e. not total flux) of the plane- 
tary nightside emission divided by that of the star, plotted as a 
function of wavelength. Using the two measurements at 8.0fj,m 
and 24-Ofj.m (black dots with error bars), we fit a blackbody curve 
through the points (gray lines) with a nightside planetary tem- 
perature ofT= 1148 ± 32ft". Open circles represent the integrated 
blackbody function across the instrument bandpasses. Filled cir- 
cles represent the same but extrapolated to the other IRAC wave- 
lengths, which allow us to conclude the nightside effect will be 
much less at lower wavelengths. 



Where we additionally define T» jCa i c as the tempera- 
ture used in the integra tion and T* b s as being equal to 
the value determined bv lTorres et al.l ((2009). Note that we 
do not fit for T„ but do allow the value to float in order 
to correctly estimate the uncertainity of Tp, n ight- The final 
analysis reveals a best-fit planetary nightside temperature 
of Tp^ight = 1148 ± 32K with x 2 = 1.12, suggesting the 
blackbody model gives a satisfactory fit for these two mea- 
surements. We note that not accounting for the response 
function of the instruments would yield a erroneous result 
of T P , night = 1120K. 

Using our derived planetary temperature, we may now 
use the blackbody function to extrapolate 72. n i g ht to other 
wavelengths and thus the nightside corrected transit depths 
for 3.6/xm, 4.5(j, and 5.8/im. The geometric transit depth will 
be given by equation (14) and the results of our analysis our 
summarized in table 3 and figure 4. 



do 



d„ 



1 — Tonight d a 



(14) 



The maximal deviation occurs for 8.0/im and is less than 
1-a for all other wavelengths. Consequently, the deduction 
of which molecules are evident from the spectrum of HD 
189733b will not significantly affected by the nightside ef- 
fect, but derived abundances will change. In future work, 
we will re-process and re-interpret the HD 189733b spectrum 
including near-infrared measurements and the nightside ef- 
fect, but a detailed re-processing of these data is outside of 
the scope of this theoretical paper. 

Combining the corrected geometric transit depths of 
table 3 with secondary eclipse measurements from the ex- 
oplanet literature, we may estimate the dayside bright- 
ness temperature in a similar manner, for completion. We 
use the secondary eclipse depth measurements quoted by 



Figure 5. Flux per unit area (i.e. not total flux) of the planetary 
dayside emission divided by that of the star, plotted as a func- 
tion of wavelength. Using the available measurements (black dots 
with error bars), we fit a blackbody curve through the points (gray 
lines) with a dayside planetary temperature of T = 1490 ± 25if. 
Open circles represent the integrated blackbody function across 
the instrument bandpasses. The nightside and dayside tempera- 
tures suggest a redistribution factor of 0. 53 



ICharbonneau et all (|2008T ) for 3.6/im, 4.5/xm., 5.8/im, 8.0/im 
and 24/xm. Note, we we do not use the 16fim data point since 
no primary transit exists for this wavelength which we can 
use to estimate to 1Zd&y A \ 2 minimization between 1000K 
and 2000K in IK steps, integrating over each bandpass in 
every simulation, reveals a best-fit dayside temperature of 
7p,day = (1490±25)K with y/x^lred = 1-45 (for comparison 
Tp,ni g ht = (1148 ± 32)K). In figure 5, we plot Tp^y/T*, in 
an analagous way as we did for the nightside in figure 4. 

Assuming the internal heating of the planet is neg- 
ligible and the planet is tidally locked, we may combine 
the two temperature estimates and employ equation (9) to 
obtain the redistribut ion factor of heat to the nightside, 
/night = 0.521 ± 0.050. iBurrows et ail (|2008l ) predict the re- 
distribution should be in the range ~ 0.1 to ~ 0.4 meaning 
our estimate may be difficult to explain theoretically. 



8 CONCLUSIONS 

We have shown that emission from the nightside of an ex- 
oplanet acts as a self-blend causing transit depths to be 
systematically erroneous, more specifically, underestimated. 
This effect is largest for infrared measurements of hot- 
Jupiter systems and can cause changes in the transit depth 
at the level of 10~ 4 , which is around the 1-a level per transit 
for Spitzer measurements of bright stars. Crucially, this is 
also the magnitude of the molecular features induced in the 
atmospheric transmission spectrum. The multiple observed 
Spitzer transits binned together raises the significance of the 
effect to 2.6-<T and for JWST we estimate the effect will be 
5-10-cr significant per transit. 

The effect may be compensated for in an accurate man- 
ner by using observations of phase curves of exoplanets to 
obtain the nightside emission flux directly and then using 
this value to correct transit depths. Naturally, this requires 
observations of both the transit event and the phase curve 
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Table 3. Using a fitted blackbody function for the nightside emission of HD 189733b, we calculate the nightside corrections to Spitzcr 



comes from 


Knutson et al\ 


\Desert et al. 





Channel 


Observed depth d bsi% 


^•night 


Corrected depth d geo ,% 


d gco — d ohs 
"d 


Measured 


8.0/^m 


2.3824 ± 0.0060 


0.1076 ±0.0020 


2.3884 ±0.0061 


1.0 


24.0/im 


2.398 ± 0.019 


0.1818 ±0.0105 


2.4085 ± 0.020 


0.5 


Extrapolated 


3.6fim 


2.356 ± 0.019 


0.0365t 


2.358 ±0.019t 


< o.i 


4.5^tm 


2.424 ±0.010 


0.0570t 


2.427 ±0.010+ 


^ 0.4 


5.8^tm 


2.436 ± 0.020 


0.0800+ 


2.441 ± 0.020 t 


$ 0.25 



in question to be at the same wavelength and ideally the 
same epoch to avoid problems associated with temporal vari- 
ability. A more practical approach may be to obtain phase 
curves at several wavelengths and use a model-fitted tem- 
perature profile to interpolate/extrapolate the required cor- 
rection for other observations. 

The effect has a strong spectral dependency and there- 
fore should be incorporated in transit lightcurve fits before 
modelling the transmission spectrum, in order to obtain the 
most accurate interpretation possible. This is particularly 
true when looking for absorption features in the infrared 
at the level of 10 -4 , which is common for many molecu- 
lar specie s already detected (for example water vapour and 
methane. ISwain et al.l (j2008|)). 
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